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Abstract 



Infrared safe differential cross sections, such, as event sliape distributions, can be 
measured over wide kinematic ranges, from regions where fixed order calculations 
, are adequate to regions where nonperturbative dynamics dominate. Such observ- 

ed ' ables provide an ideal laboratory for the study of the transition between weak and 

strong coupling in quantum field theory. This talk begins with some of the fun- 
' damentals of the perturbative description of QCD and the basis of resummation 

Q-i, techniques, followed by a brief discussion of selected topics from recent fixed-order 

^ ^ and resummed calculations. It focuses on how resummed perturbation theory has 

1^ , been used to deduce the structure of nonperturbative corrections, and to provide 

1-^^ ' a framework with which to address the transition from short- to long-distance 

dynamics in QCD. 

• 1— I 

X. 

^ . 1 Introduction 

Quantum chromodynamics (QCD) serves both as tool and background for new physics 
searches [1]. Control over QCD corrections is crucial for the interpretation of any 
experiment in which hadronic matter appears in the initial or final states, or even (as 
in the measurement of the muon's magnetic moment) only in virtual states. In addition, 
QCD is fascinating in its own right as a quantum field theory. It combines strong and 
weak coupling dynamics in any high energy experiment that is sensitive to momentum 
transfers much beyond the scale of the strong coupling, Aqcd, and every experiment 
that probes higher energies encodes the transition between partonic degrees of freedom 
at short distances to hadronic degrees of freedom at long. 

This talk approaches these issues from the partonic, short distance side, through 
the study of perturbation theory. I will argue that perturbative QCD has much to 
say about the structure of the larger theory that contains it. The discussion begins 
with a sketch of basic concepts of perturbative QCD, followed by a review of some 
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important recent developments at finite order. The extension of perturbative methods 
to multiscale observables requires resummation, and I will describe how resummed 
perturbation theory encodes nonperturbative information, with the operator product 
expansion as the classic example. Analogous reasoning is extended to jet-related cross 
sections, and examples are given of the interesting results that emerge from such an 
analysis. 

2 Some Concepts in Perturbative QCD 

2.1 Asymptotic freedom and infrared safety 

The characteristic feature of QCD is its asymptotic freedom, in which the coupling 
becomes weaker when it is measured at shorter distances. This remains a purely ab- 
stract attribute, however, until we identify observable quantities that can be expanded 
in the coupling with finite coefficients. Such quantities are said to be infrared safe [2]. 
A few cross sections, such as the total and suitably-defined jet cross sections in e'''e~ 
annihilation, are infrared safe, and may be expanded as 

n 

= ^c„(l)a^(Q) + o(i^), (1) 

where p is some integer and /x is a renormalization scale. In the second expression we 
have used the renormalization group invariance of physical cross sections. 

2.2 Factorization 

Cross sections for processes with hadrons in the initial state are not by themselves 
infrared safe. Nevertheless, if they involve heavy particle production, or more generally 
involve a large momentum transfer, Q, then short- and long-distance components of 
the cross section can be factorized [3] , usually into a convolution form, 

Q^o-phys(Q,"^) = t^SD(Q/M,as(/w)) ® fhD{l^,m) + (^-^^ . (2) 

In this expression, represents a convolution, most familiarly in terms of partonic mo- 
mentum fractions, and is a factorization scale, separating the short-distance dynamics 
of the hard subprocess, in losb, from the long-distance dynamics in /ld- /ld depends 
on various infrared, often nonperturbative scales, denoted m. For inclusive Icpton- 
hadron deep-inelastic scattering (DIS), in which Q is the momentum transfer, and for 
various other inclusive and semi-inclusive hadronic hard-scattering cross sections, /ld 
represents products of parton distribution functions. For single-particle inclusive cross 
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sections, it includes fragmentation functions. Generally speaking, 'new physics' is part 
of the short distance functions oosd, while the long-distance /ld's are universal among 
processes, making predictions possible on the basis of hypothetical couplings to new 
particles or rare processes at short distances. Factorizations also apply to many scat- 
tering amplitudes and forms factors, of which those involving the decay of B mesons 
have attracted interest recently [4]. 

Eq. (2) applies not only to factorized cross sections with hadrons in the initial state, 
but also to many infrared-safe jet cross sections, for which the scale m is large enough 
to be perturbative, but is still much smaller than the highest scale Q: Q ^ Aqcd- 
An example is the limit of 'narrow jets' in e'^e~ annihilation, where m may be taken 
as the mass of the jet [5]. 



2.3 Factorization proofs and resummation 

The proofs of relations like Eq. (2) are all based on the observation that short-distance 
processes, organized within the function wsD, are quantum-mechanically incoherent 
with dynamics at much longer length scales. This is the same observation that makes 
possible the operator product expansion. Extensions of Eq. (2) to jet cross sections and 
various decay and clastic amplitudes depend on the further observation that there is a 
mutual incoherence in the dynamics of particles receding from each other at the speed 
of light. Often, these factorization relations hold to all orders in perturbation theory 
and to all powers of \n{iJ,/Q) [6]. Interest in these issues has reemerged in analyses of 
high energy processes in the language of effective theories [7]. 

Whenever there is factorization, there is some form of evolution and resummation. 
To see why, we observe that a physical cross section cannot depend on the factorization 
scale, 

= In o-phys(<5, rn) . (3) 

From this it follows that the variations of long-distance and short-distance functions 
with respect to /x must compensate each other, through a function of the variables that 
they hold in common: 

d d 
'^djj. ^^•^ld(m,"^) = -Pi^sil^)) = IntusD {Q/li,as{lj)) ■ (4) 

In deep-inelastic scattering, and other inclusive hard-scattering cross sections, where 
wsD and /ld are connected by a convolution in momentum fractions, the functions 
P{z,as) are the DGLAP splitting functions [8]. 

By the same token, wherever there is evolution, there is a resummation of per- 
turbation theory, generated by solving the evolution equation. Again thinking of DIS, 
we may fix /x = Q in Eq. (2), compute wsd as a power series in as{Q), and derive the 
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Q-dependence of the cross section by solving (4) for the functions fLD{Q,n^)- 

lncrphys(Q,"i) ~ UJsD{l,as{Q)) 

®exp||^'^^P(a,(M'))| fLD{q,m). (5) 

The most impressive successes of this approach are in observables with a single hard 
scale. Inclusive DIS is the exemplary case, as illustrated in Fig. 1 for the structure 
function F2{x,Q), which is proportional to the cross section to produce an hadronic 
final state of total invariant mass squared equal to Q^{1 — x)/x, at momentum transfer 
Q. 
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Figure 1: Data for F2{x,Q) with QCD fit based on evolution, from Ref. [9]. 

As mentioned above, however, only a limited number of observables depend on a 
single hard scale. A more typical example is the distribution of p = (mj/Q)'^, from 
e'''e~ annihilation dijet events, with mj the mass of the heavier jet in the center-of-mass 
(cm.). The experimental distribution in p at the Z pole is illustrated in Fig. 2, taken 
from [10]. In this case, the two scales are the cm. energy, Q = y/s and the mass of 
the 'heavy' jet. The data is shown as a function of p = m'j/Q'^, and even resummed 
perturbation theory (next-to-leading logarithm in rrij/Q) is not adequate to follow it. 
It is necessary to include power corrections in the lighter scale: ^/mFj (as organized in 
a shape function, see below), which become more and more important in the 'exclusive 
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limit' of light jet masses. Notice that most events are in this range; going one step 
beyond the fully inclusive cross section, we are required to find a description with 
nonperturbative as well as perturbative input, even though the observable in question 
is infrared safe. As we decrease mj, nonperturbative effects vary continuously from 
a few to tens of percent. Thus, within this single set of data, at a single energy, we 
can select events whose formation is sensitive to an adjustable mix of perturbative and 
nonperturbative dynamics. This is an ideal testing-ground for quantum field theory. 



Figure 2: Heavy jet distribution at the Z pole, from Ref. [10]. The dashed line is 
resummed perturbation theory, and the solid line a fit with a shape function. 

I shall discuss below recent attempts to organize nonperturbative corrections to 
infrared safe cross sections starting with perturbation theory, in a manner similar to 
factorization theorems like Eq. (2). To do so with confidence, we should attempt 
to organize perturbation theory to all orders in the coupling, the level at which we 
understand Eq. (2). First, however, let us touch on some of the very substantial progress 
of the past few years that is based on finite-order calculations. 

3 Toward a Two-Loop Phenomenology 

Fig. 3 is a cartoon of the progress of fixed-order calculations in perturbative QCD over 
the past three decades, with tree amplitudes under control around 1980, one-loop am- 
plitudes around 1990, and two-amplitudes since 2000 (early two-loop advances were 
reported [11], with recent reviews in [12]). In the years following their calculation, am- 
plitudes at the tree and one-loop level were implemented into factorized cross sections 
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like (2). With the hard-scattering computed at tree level, the cross section is termed 
'leading order' (LO). At one loop, it is 'next-to-leading order' (NLO), and so on. Now 
we are at the dawn of an age of NNLO. It is still the dawn, because loop diagrams suffer 
from various infrared divergences, which must be controlled by adding gluon emission 
diagrams and constructing useable, factorized hard-scattering functions. It is a big step 
from a set of (dimensionally) regulated amplitudes to a factorized cross section (2) that 
can be compared to experiment. Still, with the two-loop amplitudes in hand, the goal 
of NNLO phenomenology is certainly attainable. 




Figure 3: Progress in low-order calculation. 

Why go to the trouble of these very challenging calculations? With as{Q) ~ 0.1 
at Q ~ 100 GeV, computations to order can reach the accuracy of one percent, an 
admirable goal in its own right for the strong interactions. At the same time, accuracies 
of a few percent will enable a new class of investigations of QCD, by covering kinematic 
regions where nonperturbative power corrections vary from a few percent to tens of 
percent, as we have seen above. Beyond this, one percent accuracy will make possible a 
new generation of determinations of the strong coupling itself, the uncertainty in whose 
precise value is a major limitation on our ability to extrapolate the standard model into 
the realm of new physics [13]. Recent mileposts on this road include the calculation of 
0{as^) 3-jet amplitudes in electron-positron annihilation [14]. 

For factorizable jet and other cross sections at n loops in the hard scattering, it 
is desirable to have n + 1 loops in the splitting function P{z, as), in Eq. (4). Progress 
toward P{z, Us) to has included numerical estimates [15], leading to studies of deep- 
inelastic scattering at this level [16], as well as a growing development of the methods 
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necessary for full analytic expressions [17] at three loops. The coming decade should 
see the full realization of the second step of Fig. 3. 

The long-distance components of factorized cross sections are parton distributions 
(PDFs) and fragmentation functions. Unpolarizcd parton distributions implemented 
with NLO calculations of hard scattering have reached a stage of refinement over large 
ranges of parton fraction x. These refinements are of great significance, because they 
are a limiting factor in our ability to extrapolate known cross sections to the scales at 
which new physics signals may appear. This said, there have still been 'surprises', in 
just the past few years, such as the unexpected behavior of sea antiquarks as measured 
by experiment E866 at Fermilab [18], and promptly incorporated into the following 
round of parton distribution functions [19, 20]. 

Perhaps the most important recent development in the technology of PDFs is 
the effort to quantify their uncertainties [21]. This is a daunting task, because any 
set of PDFs for all of the quarks, antiquarks and the gluon requires the specification 
of a substantial set of parameters. Generally, in a 'global' fit, these parameters are 
fixed by comparison with a set of existing data. Any estimate of the uncertainty in a 
calculated process depends in a unique way on all of these parameters, each of which 
is sensitive to both statistical and systematic errors in the original data set. What 
we are interested in is not so much the uncertainties in the PDFs themselves, but 
what these uncertainties mean for cross sections that we would like to calculate. To 
determine this, it is necessary to specify, not just a 'best fit' for the PDFs, but a set of 
degrees of freedom accompanied by correlated error estimates. With this information, 
it is possible to create an ensemble of PDFs, selected according to their statistical 
likelihood. Each element of such an ensemble provides an independent calculation 
of a factorized cross section (2), and the distribution of these results determines the 
desired uncertainty. Needless to say, this procedure is computationally-intensive in all 
of its variants, depending on the choices of parameterization, the analysis of parameter 
uncertainty, and the implementation of the ensemble-creation process. The rewards, 
however, are substantial, and the effort is mandatory for extrapolation to the range of 
the highest-energy hadron accelerators. 

4 Resummation 

Fixed-order calculations work best for single-scale observables, with the jet cross sec- 
tions in hadronic collisions offering a fine example, both of the range of applicability, 
and the uncertainties involved [22]. The cross section for the production of bottom 
quarks is an example involving more controversy, and possible recent progress [23]. 
Any finite-energy calculation, however, leaves the nagging question of what happens at 
yet higher order. For single-scale inclusive cross sections, the comparison of LO with 
NLO and NNLO, when available, can give confidence or shake it. (See, for example, 
recent NNLO calculations of inclusive Higgs production at hadron colliders [24].) Nev- 
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ertheless, there are certain general features of perturbative corrections that demand 
study at all orders, even in one-scale cross sections, and in most cases two-scale cross 
sections cannot be controlled at all without a classification and summation of such 
corrections at all orders. Examples of these two possibilities are threshold and Qt 
resummation. 



4.1 Threshold resummation 



Threshold resummation addresses a specific set of corrections at all orders in pertur- 
bation theory for factorized cross sections. Consider a generic factorized cross section 
for the inclusive production of an observable final-state object F of invariant mass Q, 
written in terms of partonic fractions as 



Q^UAB^FiQ) = fa/A{Xa) ® fb/Bi^b) ^ ^ab-^F+X Z 



(6) 



with the fc/H PDFs for parton c in hadron H. The hard-scattering cross section 
depends on variable z, which is the ratio of to the squared invariant mass of the 
pair of partons that initiate the process ab F + X, where X denotes additional 
radiation in the inclusive final state. Of special interest is the limit z — 1, which is 
called partonic threshold. At partonic threshold, the partons a and b have only enough 
energy to produce F{Q), with nothing left over for radiation. This is a special limit 
in field theory, because a hard scattering favors processes with copious radiation. The 
consequence of this mismatch is that z = 1 is associated with singular but integrable 
distributions, such as 



(r) 

ab^F 



1 

H 



In' 



2r-l 



{1-z) 



(7) 



+ 



where a plus distribution [5(0;)]+ with g[x) singular at x 
function) by its integrals with smooth functions: 



1 is defined (like a delta 



['dx[g{x)]+f{x)= [' dxg{x)[f{x)-f{l)]-f{l) j' dxg{a 

Jy Jy Jo 



(8) 



In Eq. (6), the singular distributions are in ujab^F+x{z), which has the information on 
perturbative radiation, and the smooth functions are the PDFs. While all of the plus 
distributions are integrable, they become more and more singular at higher orders, and 
it is natural to try to estimate their overall effect beyond the lowest orders. 

We are able to resum singular plus distributions because at partonic threshold the 
hard-scattering process returns to the elastic limit that we have referred to above, in 
which final state radiation is grouped into well-collimated jets, associated either with 
the fragments of the incoming hadrons, or with a set of nearly light-like outgoing jets. 
Such a momentum configuration is illustrated by Fig. 4 for the production of two jets. 
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Figure 4: Elastic kinematics and refactorization scale at partonic threshold. 



At partonic threshold, we may implement the sequence of reasoning that we de- 
scribed in Sec. 2. Quanta are part of well-collimated jets, moving relatively at the speed 
of light, or are very soft. We may (re)factorize the partonic cross section in this limit, 
as indicated in Fig. 4 by the circle and 'refactorization scale', which leads to new 
evolution equations, and hence to resummation. For the hard scattering of Eq. (6), we 
find expressions whose expansion generates the terms identified in Eq. (7) to all orders 
in perturbation theory [25]. The resummation is most conveniently exhibited in the 
space of Mellin moments of Uab^F+Xi where they exponentiate as 



dzz^ ^UJab^F+x{z) 



exp 



exp 



i=a,b 



1 _ 1 



E / dz- 



Z J{l-z)2Q 



2 m 



E 

i=a,b 



drr? 



Ai{as{m)) In 



fNm\ 

\q) 



(9) 



in terms of anomalous dimensions Ai{a^ = Ci{as/'n') + . . ., with = Cp for i = q,q 
and Ca for gluons. Eq. (9) generates the leading logarithms at each order; the complete 
expression involves process-dependent nonleading logs. At each order in the exponent, 
one of the logarithms can be evaluated explicitly, as shown, but in the remaining integral 
the strong coupling runs down to arbitrarily low scales. We will return to this intriguing 
feature of the resummed cross section below. There are various ways of dealing with 
this low-scale limit [26, 27, 28, 29], and although the issue is far from settled, the 
bottom line seems to be that the effect of singular distributions at partonic threshold 
is a modest enhancement, along with a reduction of dependence on the factorization 
scale. 
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4.2 Color flow in resummation and at fixed order 



The expression (9) involves radiation that is nearly collinear to the incoming and out- 
going partons at the hard scattering, and is completely incoherent among these sets 
of lines. It organizes the leading logarithms in N, and a portion of subleading correc- 
tions. Beyond leading logarithm, however, the active partons also radiate soft gluons 
coherently, in a manner that depends on the exchange of color at short distances - 
schematically, inside the dashed circle labelled /i' in Fig. 4. At this level, thresh- 
old resummation has an interesting interplay with the one- and two-loop calculations 
[11, 30, 31] represented by Fig. 3. 

In the resummed cross section, the effects associated with coherent soft-gluon 
emission are single- logarithmic in the moment space of (9), 



exp 



dm 

Q/N m 



(10) 



where the A's are eigenvalues of a set of anomalous dimension matrices, which can be 
determined by separation of variables as in Eq. (4). They are matrices because the 
short-distance function of Fig. 4 (inside the circle) can be thought of as a vector in the 
space of color exchange. For example, when the active partons are all gluons, there 
are as many as eight independent color exchanges at short distances. As fi' changes, 
gluons that are exchanged near scale ji' enter or leave the hard scattering. These effects 
are reflected in the NLO amplitudes of Fig. 3, and influence the infrared structure of 
NNLO in an essential way [30, 31]. 

A striking example is the gg gg anomalous dimension matrix, which controls 
single logarithms in threshold resummation for gluonic jet production [32], with A^^, the 
number of colors: 
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(11) 

1 — cosd*), U = ln^(l -|- cos^*) in terms of the cm. scattering angle 
6*. The same matrix appears in explicit two-loop calculations of elastic gluon-gluon 
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amplitudes [11], where it controls the single poles in dimensional regularization. The 
all-orders nature of this relationship has been demonstrated in [31]. 



4.3 Qt and joint resummation 

If threshold resummation strengthens our confidence in cross sections that we can 
calculate to fixed order, there is a class of hadronic cross sections for which resummation 
is absolutely necessary even to get started [33, 34]. The most familiar of these are the 
hadronic electroweak annihilation cross sections at measured transverse momentum 
{Qt) for the production of an electroweak boson B{Q) (virtual photon, W, Z or Higgs, 
of mass Q). In the underlying processes qq B{Q) or gg B{Q) (the latter through 
a fermion loop), the limit Qt ^ is singular in much the same way as 1 — z ^ at 
partonic threshold, Eq. (7), with plus distributions in Qj, up to (l/Qj^) ln'^^~^{Qj./Q'^) 
at order ctg'^. Indeed, these two limits are closely connected and involve the same 
underlying factorization, and the resummation necessary for logarithms in transverse 
momentum is compatible with threshold resummation [28, 35]. Again, the result for 
the factorized/resummed cross section exponentiates in transform space, in this case a 
combination of Mellin (N) and Fourier (6, 'impact parameter') transforms. The 'jointly 
resummed' physical cross section is then approximated by a double inverse transform, 

dfT'"^' _ ^aa^B{Q) f dN f d^b ^^Q^,^ 



dQ^dQ^T a Jc 27rz Jc, {lirf 

X C„(Q, h, N) e^°^ C-a{Q, b, N) , (12) 

where the PDFs are included in the factors 

Ca{Q,b,N,n) = Y.Ca/j{N,as{fi)) fj{N,Q/{N + bQ)) . (13) 
i 

The leading logarithmic corrections are of much the same form as in threshold resum- 
mation, Eq. (9), 

EaaiN,b,Q) ~ 2 /'^ ^AMm)) In(^) , (14) 

J[Q/{N+bQ)]^ \QJ 

in terms of the same anomalous dimensions Aa{as)- A resummed cross section cal- 
culated from (12) is compared to Tevatron data for Z production in Fig. 5. Similar 
good fits are found on the basis of Qr-i'esummation alone [33, 34]. We note in the 
figure a general resemblance to the heavy jet mass distribution of Fig. 2, where now 
the exclusive limit is at Qt = 0, which has the kinematics of quark-antiquark annihila- 
tion without gluon radiation. There is a need for modest [34, 36] but nonzero 'power' 
corrections, as shown in the figure. 
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Figure 5: Jointly-resummed cross section da/dQT for Z production at the Tevatron, 
from Ref. [36]. 

5 Power Corrections and Event Shape Functions 

Examining Figs. 2 and 5 we see that for differential distributions, power corrections 
are necessary to bring even resummed cross sections into agreement with data in the 
exclusive limit. In the range where the kinematics are far from elastic, fixed-order 
perturbation theory is enough, but as we climb the curves of these figures toward the 
elastic limit, the cross section rises and we move into a region where first resummation, 
and then nonperturbative corrections dominate. The theory itself knows no strict 
boundary between perturbative and nonperturbative regimes. It is therefore natural 
to analyze perturbation theory to see if it gives signals of its own incompleteness, 
and to seek hints on how to supplement perturbative predictions with nonperturbative 
information. The ultimate goal is to construct a single theory that bridges the gap 
between short- and long-distance dynamics. 

5.1 The operator product expansion 

The approach that I will describe below begins with the simplest one-scale problem, 
the total cross section in e'''e~ annihilation [37]. The total cross section is related by 
the optical theorem to the imaginary part of the vacuum polarization for electroweak 
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currents, 



X e-"i-^Q\Tf{Q)j^{x)\Q), 



(15) 



whose behavior at large = is determined by an operator product expansion: 



Here, I have shown only the contributions of the leading (identity) operator and of the 
gluon condensate F'^. Neglecting masses, perturbation theory alone would give only 
the leading term, even though the full theory demands both (and more). The question 
we ask is whether we could have discovered the need for the nonleading operators from 
perturbative QCD alone. 

The answer is yes, and the reasoning is illustrated in Fig. 6, which shows graphical 
configurations where loop momenta k are vanishingly small for the vacuum polarization 
and for the gluon condensate treated as a local operator. Diagram-by-diagram, each 
such configuration is a tiny 'corner' of loop momentum space, where one or more 
loops has low momentum, giving a finite contribution to the infrared safe integrals. 
For example, if the loop momentum is required to have k"^ < k^, for some cutoff k, 
this contribution is suppressed by {k,'^ /Q'^) as{fi) compared to the full diagram. It is 
thus indeed power suppressed. Using gauge invariance and the renormalization group, 
however, we can show that summing over the self-energies indicated in the figure has 
the effect of making the coupling run and of giving an even greater suppression. For 
the vacuum polarization, we find momentum-space integrals of the general form 



with H(Q) an infrared safe function free of all contributions where a single loop mo- 
mentum vanishes. H{Q) represents the lower, contracted part of the diagram on the 
left of Fig. 6. Up to this factor, the expression for the right-hand figure, corrections 
to the condensate, is identical. A generalization of Eq. (17) is valid to all orders in 
perturbation theory, and controls all logarithms in the loop momentum k, up to the 
next power, Q^^ [37]. We will call an expression at this level of accuracy, all logs at a 
given nonleading power, an 'internal resummation' at that power. 

From (17) we learn first, that the internally-resummed contributions are singular 
for very small fe^, from high order diagrams. Correspondingly, it is easy to verify that 




(16) 





(17) 
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the nth order of the reexpansion of Eq. (17) in as{Q) is proportional to n\. We also learn 
that these infrared-sensitive contributions in the vacuum polarization generate exactly 
the perturbative expansion of the gluon condensate. We conclude that the perturbative 
vacuum polarization requires the presence of a nonperturbative gluon condensate, which 
should replace its misleading high-order infrared behavior. 




Figure 6: Sources of infrared sensitivity in perturbation theory for Co (left) and (-F^) 
(right). 

The foregoing steps may be summarized as an 'axiom of substitution', in which 
the infrared-sensitive, nonconvergent portion of an internally-resummed perturbative 
quantity is replaced by a matrix element with the same perturbative infrared expansion: 

CpT ^ C'p'^{Q^/fi\K/Q,as{fi)) 

+ -^CMQ,^) c^s{0\F\0)\0)iK). (18) 

Here, Cp*^ is perturbation theory subtracted in its infrared-sensitive limit. As antici- 
pated, the nonconvergence of perturbation theory implies the need for a new, infrared 
regularization. For the vacuum polarization the cost is a new nonperturbative param- 
eter {{F"^)), necessary to define the theory at the level of l/Q^ relative to the leading 
behavior. But, because it is implicit in perturbation theory, this new parameter is at 
the same time a reward of the analysis. 

5.2 Power corrections in semi-inclusive cross sections 

The viewpoint described above for the operator product expansion has been imple- 
mented, in various ways, to a variety of infrared safe cross sections and hard-scattering 
functions [38]. Of particular interest are resummed hard scattering functions like Eqs. 
(9) and (14), which exhibit the same integral over the scale of the coupling as for the 
vacuum polarization (17). The application of these methods to event shape cross sec- 
tions [39] sheds light on the transition from short- to long-distance dynamics in the 
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formation of final states in QCD. This section will sketch the ideas behind these ap- 
plications, and also touch on internal resummation in this context. Such questions are 
most easily addressed for event shapes of two-jet cross sections in e^e~ annihilation, 
although their application is much wider [40]. 

We consider inclusive cross sections defined by event shapes, S{{pi}), 

^ = ^J2[ \Mn{{Pr})\'6 {S{{p.}) - S) , (19) 

at cm. energy Q, computed from QCD matrix elements M„({pi}) integrated over n- 
particle phase space, PS{n). This event shape cross section is infrared safe if S{{pi}) is 
a smooth function of final state momenta, pi , that is unchanged by soft parton emission 
and collinear rearrangements, as illustrated in Fig. 7: 

S{. ..pi . ..pj-i,api,pj+i ...) = S{...{1 + a)pi . . . . .) . (20) 

This condition makes it possible to compute the cross section perturbatively, because it 
deemphasizes the long-time evolution of the system, over which the nearly degenerate 
states of Fig. 7 mix. 

The study of event shapes in QCD, and more generally energy flow [41], continues a 
long tradition in field theory. We may draw a not too far-fetched analogy between how, 
late in the nineteenth century Poynting discovered the rules for energy and momentum 
flow in classical electrodynamics by considering the slow discharge of a condenser [42], 
and how, late in the twentieth century, the short-distance structure of QCD was re- 
vealed in the rapid neutralization of a quark-antiquark color dipole produced in e"'"e~ 
annihilation [43]. 




ap 



Figure 7: Soft emission and collinear rearrangements. 

We restrict ourselves here to a specific set of weights, appropriate to dijet events 
in annihilation cross sections [44], 

'5-W = ^ E ^.T e-l^'Ki-") , (21) 

^ all iGn 
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where the transverse momentum kiT and pseudorapidity rji of each particle momentum 
ki is defined relative to the axis that minimizes the specific function 5o(n) for state 
n. The variable a takes on any value less than two. The particular choice a = is 
Sq = 1 — T with T the thrust, while o, = 1 is known as the jet broadening [45]. 

For these dijet event shapes, sensitivity to long times can be adjusted by studying 
the Laplace transform [39, 45], 

aa{N) = ^ f'dSae-''^^^. (22) 

C^tot JO a^a 

Choosing N large, we select states with small Sa, which are sensitive to long-time 
behavior (roughly, times of order 1/SaQ), corresponding to low jet masses in dijet 
events. As illustrated by Fig. 2, however, power corrections are important at low jet 
masses. We thus ask what the reasoning described above for the operator product 
expansion can tell us about the transformed cross sections, especially at large A'' in Eq. 
(22). 

At low values of Sa, these differential event shape cross sections show the same 

double logarithmic behavior discussed above for threshold rcsummation, and obey a 
similar resummation in Laplace transform space [44, 45], which is given to leading 
logarithm, accompanied by effects of the running coupling, by 

{rQ'^ /7t)2 fPr/Q rill 

Jo pt Jp^t/q^ y 



(23) 



in terms of the same Sudakov anomalous dimension as above, Aq{as) = Cpiois/T^) + . . .. 

We can use Eq. (23) to suggest the structure of power corrections in differential 
cross sections, as described above for the total annihilation cross section. To do so, we 
split up the Pt integral into 'perturbative' and 'soft' ranges, separated at a factorization 
scale K [46] . If we expand the exponent of (23) in the soft range, we find 



£^(-f)7"'#pM(«.w) 

I- a nn\ \ QJ Jo Pt 



1 



'Q 



= lnaa,PT(iV, Q, k) + In f\{N/Q) + O ( ] . (24) 



N 



In the second line we organize the entire contribution from the soft region into a single 
function fsai^/Q)^ with corrections relatively suppressed by powers in (assum- 
ing for now that a < 1). Taking the inverse transform of (24) we arrive at an expression 
for the physical cross section as a convolution of the perturbative cross section with 
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what is known as an event shape function [46], 

do a Jo do a 

The event shape function fg^ is independent of Q, so that a fit to data at Q = mz is 
sufficient to predict the differential cross section at any Q, including all integer powers 
of N/Q in moment space [10], which translates to all integer powers of 1/SaQ at fixed 
Q. This is the kind of formalism necessary to confront data like those in Fig. 2 above. 
Analogous functions find important applications in the physics of b quarks [47]. A 
recent example appled to the heavy jet mass, closely related to a = 0, is shown in Fig. 
8 from Ref. [48]. Beyond the a = 0, the shape functions derived in this fashion for 
choices of a, 6 < 1 possess a surprisingly simple relation, 

In fsAN/Q) = ^ In fs, {N/Q) , (26) 

with corrections as in (25). For 2 > a > 1, on the other hand, we encounter the 
interesting situation noted first in Ref. [49], where the leading power correction becomes 
fractional, and larger than l/Q. 




0.1 0.2 0.3 p„ 



Figure 8: Event shape-based fit to the heavy jet mass at various energies, from Ref. 
[48]. 

The next question we should ask, however, is whether Eqs. (23) and (26) hold 
beyond leading and (suitably generalized [45]) next-to-leading logarithmic resumma- 
tion, and in particular whether these event shapes can be internally resummed to all 



SaQ 



2-a 



(25) 
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orders in 1/SaQ- Rather than try to answer this question for the full theory, we will 
follow Ref. [46], and consider the 'eikonal' cross section, in which the primary quark- 
antiquark pair and the partons collinear to them are all replaced by nonabelian phase 
operators. These are Wilson lines, represented by products of exponentials of the gauge 
field ordered along lightlike paths, 

T^ = r(^$(f(cx),0)cl.g(oo,0)j, (27) 

with 

$g (A, x)=V exp (^-ig £ dX'Pi ■ A^'^ (A'A + x)^ , (28) 

where yl^*-* is the gluon field in representation i. In these terms, the eikonal cross 
sections at fixed values of event shapes Sa may be written as a sum over states, n: 

1 ^^eikonal 

= J2S{Sa- Sain)) {0\W^ In) (n| W |0) . (29) 

Ctot (l^a n 

These cross sections have the advantage that many of their properties may be analyzed 
to all orders in perturbation theory using combinatoric methods [50] . Here I will simply 
quote the result of such an analysis. 

As in Eq. (23), the all-orders eikonal approximation to the cross section can be 
written in transform space in terms of a function in which runs with an integration 
variable, which we may think of as a transverse momentum, 

)2 j2 



In a„(iV, Q)= f"^ ^ As^ (pT/fi, Q///, «.(/x), N) , (30) 
Jo Pt 



where the function Asa is infrared safe and renormalization group invariant {dAsa /dfJ- = 
0). As an example, consider the case a = 1, jet broadening, for which we can write 



-^51 {pT/fJ',Q/fJ',Ols in) , N) 

lps{n) ( ' • h 



n 



(31) 



where the ki denote particle momenta in final state n, and where rjn is the rapidity of the 
total momentum of state n. Boost invariance requires that the functions VV„ [50] depend 
only on the arguments shown. For other event shapes, such as the thrust, slightly more 
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elaborate expressions are necessary. In each case, however, a meaningful factorization 
of the infrared, strong-coupling dynamics from the perturbative region is possible. Eq. 
(31), although valid for all orders and logarithms, applies only to the eikonal cross 
section. Nevertheless, it appears that the way is open for the all-orders analysis of the 
relation between perturbation theory and power corrections in event shape functions 
and related parameterizations, beyond the limitations of next-to-leading logarithm. 

6 Summary 

Short-distance, single-scale observables in QCD are reasonably well understood, which 
is one reason why we have such confidence that QCD is a true theory of the strong 
interactions. The developing technology of next-to-next-to-leading order, coupled with 
a sophisticated approach to parton distribution uncertainties, should make possible 
precision for selected observables at the level of one percent. 

At the same time, we are far from a full understanding of multi-scale observables, 
which lead us from perturbative to nonperturbative dynamics. Power corrections and 
event shape functions for infrared safe differential distributions may be identified by 
an extension of the reasoning that leads to the operator product expansion for the 
total e+e~ annihilation cross section. By exploring concepts like internal resummation, 
we may find a perturbative window to the formation of final states, and gain insight 
into the transition between short- and long-distance degrees of freedom in quantum 
chromodynamics. 
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